Crystallization of hard-sphere glasses 



E.Zaccarelli, 1 C.Valeriani, 2 E.Sanz, 2 W. C. K. Poon, 2 M. E. Cates, 2 and P. N. Pusey 2 

Dipartimento di Fisica and CNR-INFM-SOFT, Universita' di Roma La Sapienza, Piazzale Aldo Mow 2, 00185, Roma, Italy 
SUPA, School of Physics and Astronomy, The University of Edinburgh, 
James Clerk Maxwell Building, The King's Buildings, Mayfield Road, Edinburgh, EH9 3JZ, UK 

(Dated: October 20, 2009) 

We study by molecular dynamics the interplay between arrest and crystallization in hard spheres. 
For state-points in the plane of volume fraction (0.54 < <f> < 0.63) and polydispersity (0 < s < 0.085), 
we delineate states that spontaneously crystallize from those that do not. For non-crystallising 
(or pre-crystallization) samples we find isodiffusivity lines consistent with an ideal glass transition 
at (j> g ~ 0.585, independent of s. Despite this, for s < 0.05, crystallization occurs at <f> > <f> g . 
This happens on timescales for which the system is ageing, and a diffusive regime in the mean 
square displacement is not reached; by those criteria, the system is a glass. Hence, contrary to a 
widespread assumption in the colloid literature, occurrence of spontaneous crystallization within a 
bulk amorphous state does not prove that this state was an ergodic fluid rather than a glass. 
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Pioneering computer simulations in the 1950s pre- 
dicted that a system of monodisperse hard spheres (HS) 
should crystallize at high enough volume fraction, <f> [lj]. 
Thirty years later, a full phase diagram (from fluid den- 
sities up to random close packing, 4>rcp ~ 0.64) was 
measured in a suspension of hard-sphere PMMA colloids 
Q. Alongside the fluid and crystal phaseSjthis reported 
the formation of a glass at <p > <p g ~ 0.58 @, The ex- 
istence and nature of the HS glass remains controversial 
@, H, H 0] : for instance, a recent report @ suggests that 
HS suspensions remain ergodic well beyond <f> = 0.58. 
Ergodicity depends on observation time, so this finding 
does not rule out colloidal glasses on timescales of several 
hours. (It does rule out an ideal glass transition, in which 
the structural relaxation time for cf> > 4> g is strictly infi- 
nite, as predicted by mode-coupling theory (MCT) [§].) 
For many years, a number of colloid physicists have held 
that, because particles in a glass cannot rearrange dif- 
fusively, crystallization within a homogeneous bulk glass 
cannot proceed unless it is seeded with pre-formed nuclei 
This view is stated explicitly in Q (but see [HI); 
in combination with the observation that hard sphere 
colloids at 4> > <f> g show rapid bulk crystallization in mi- 
crogravity, it has led Chaikin and others to assert that 
such colloids do not have a glass transition at all 12 1. 



An important factor in studies of all these phenomena 
is polydispersity. Fractional standard deviations in parti- 
cle size, s, of a few percent are inevitable experimentally. 
Some authors hold that polydispersity is essential to the 
formation of a hard sphere glass [13, [lj] because a pu- 
tative monodisperse sample would crystallize too rapidly 
[lit . (Here, as for glasses generally, the contest between 
vitrification and crystallization is usually portrayed as 
happening during rather than after a quench (lf3|.) In col- 
loids, polydispersity could play two distinct roles. First, 
it might influence cage formation, shift ing any glass tran- 
sition point. However, computational [17l.ll8l] and exper- 



imental [19( | work seems to suggest that modest polydis- 
persity (s < 0.1) has little effect on structural relaxation 
times. On the other hand, numerical and exper- 

imental [HI, HH studies have shown that crystallization 
is strongly suppressed even at small finite s. Hence a 
second possible role of polydispersity is to destabilize di- 
rectly the crystal. Indeed work on the equilibrium HS 
phase diagram [22|, 0, 0, predicts that beyond a 
'terminal' s value (~ 0.07) there is no stable crystal of 
the same composition as the fluid. Crystallization then 
requires size fractionation [2BI ], which involves transport 
over large distances. 

In this Letter we present molecular dynamics (MD) 
simulations that probe in detail the effects of polydisper- 
sity on the HS glass transition and on crystallization. We 
find almost no effect of polydispersity on the dynamics of 
diffusive rearrangement within the amorphous state near 
the glass transition, but confirm also that polydispersity 
strongly suppresses crystallization [l3|, • More strik- 
ingly, we find that for low polydispersity (s < 5%), hard 
spheres crystallization can occur spontaneously within a 
bulk amorphous state that is plainly a glass (no particle 
diffusion) on the timescale of observation (see Fig. 3 be- 
low). Note that, just as in colloids [1, [l(| non-diffusive 
pathways to crystal growth -from a pre-seeded nucleus, 
within a bulk glass- are well-established for conventional 
(microscopic) glass-formers [26j |. Experimental evidence 
for crystallization in unseeded homogeneous bulk glasses 
exists [27J, but is much less clear: surface devitrification 
and residual diffusion are both hard to rule out. Our sim- 
ulations establish that crystallization without diffusion is 
indeed possible in the simplest of all glasses: bulk homo- 
geneous hard spheres. We plan further detailed analyses 
to elucidate the mechanism; future colloidal experiments 
could equally illuminate the microphysics involved. 

Our MD study addresses N — 2000 hard spheres in an 
NVT ensemble, with 0.54 < < 0.63 and < s < 0.085, 
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using an event-driven algorithm for particles interacting 
via hard potentials. To approximate a bulk system, we 
use periodic boundary conditions in a cubic box (volume 
V). Mass, length and time are measured in units of the 
particle mass m, the mean particle diameter a, and the 
thermal timescale to = ^(m/fcsT) 1 / 2 , where fegT = 1. 
Polydispersity is represented by a 31 component discrete 
Gaussian distribution of standard deviation s; the vol- 
ume fraction obeys cf> = irNa 3 /6V. At each 'state point', 
defined by given values of <f> and s, we run simulations of 
duration tM = 10 5 ; to improve the statistics we perform 
5 independent runs. Each run is started from a different 
initial configuration; those containing ordered seeds are 
eliminated by requiring an initially low percentage (less 
than 5%) of 'crystalline particles'. (We check that any 
remaining patches of initial order do not act as nuclei.) 
Such particles are defined locally u sing a rotationally in- 
variant bond order parameter de [281 ] . with cutoffs set 
at r c — 1.4 (to define each particle's first neighbors), 
d c = 0.7 (a first criterion that identifies candidate crys- 
talline particles as having de > d c ) and N c — 6 (a second 
and more stringent criterion requiring such a particle to 
have at least N c neighbors that also satisfy d§ > d c ). For 
additional details and results of the simulations, see [29j |. 

For each state point we track the pressure, and the 
fraction of crystalline particles, over time: a state point 
is declared 'crystalline' if at least one of the five simu- 
lation runs develops a proportion of crystalline particles 
exceeding 85%, within t^.j- The resulting stability dia- 
gram in the (</>, s) plane is shown in Fig. [TJ This broadly 
agrees with, but vastly extends, the much earlier simula- 
tion results of [131 ]. 

Fig. Q] shows that crystallization takes place in a well- 
defined region of the stability diagram. The maximum 
polydispersity observed, s ~ 0.07, almost coincides with 
the predicted limit of existence of a non- fractionated solid 
[24j ]. Polydispersity is relevant on both the low and the 
high-0 borders of the crystallising region, narrowing to 
a similar extent the crystallizing range from both sides. 
The \ow-(f> part of the boundary almost tracks the pre- 
dicted density of the fluid at coexistence, with a fixed 
supersaturation A(f> w 0.05. Thus the effect of s on crys- 
tallization at low (j) can be accounted for simply by the 
shift of thermodynamic equilibrium boundaries (which 
reduces supersaturation at a given (f>) , combined with an 
assumption that a nucleation barrier is only surmount- 
able for sufficient A(f> [2(J. On the other hand, as stated 
previously, suppression of crystallization at high 4> could 
stem from slowed kinetics in the amorphous phase, or di- 
rect destabilization of the crystal. To distinguish these al- 
ternatives, we next assess the impact that polydispersity 
has on structural relaxation dynamics within the amor- 
phous state. 

For this we require quantitative data on long-time dif- 
fusion, and so must focus on those state points which re- 
main disordered during our simulation time. For each of 
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FIG. 1: Stability diagram in the ((j), s) plane, reporting state 
points that crystallize (black crosses) in at least one of five 
runs, and state points where all runs remain disordered until 
tM (red circles). We also show as triangles isodiffusivity lines 
(left to right: D = 5 x 10" 3 (blue), D = 10" 3 (violet), D = 
5 x 1CP 4 (cyan), D=5 10 -5 (orange)); to the right of these lie 
the MCT glass transition for a binary mixture of matched s 
(maroon, open triangles), and the extrapolated ideal glass line 
from our data (magenta), with a dotted continuation into the 
crystal as guide to the eye. The phase boundaries allowing [24[ 
and not allowing 25] for size fractionation are shown as solid 
(turquoise) and dashed (green) lines respectively. 



these we calculate the mean square displacement (MSD) 
and extract from its long-time limit the self-diffusion co- 
efficient D (averaged over 5 runs). In Figure [2] (top) 
we plot the run-averaged MSD curves {(r 2 (t))) for non- 
crystallizing state points at all s with 0.54 < <fi < 0.58. 
Of these, only for s = 0.07 and s — 0.085 can we study 
the full </> range directly. But for no </> is there evidence 
of significant s-dependence in the diffusive rearrangement 
kinetics. Instead we recognise at all s the same short and 
intermediate time behavior, and observe the emergence 
of the same plateau, typical of systems close to dynamic 
arrest, for <fi > 0.57. 

Fig. [U (bottom) plots the extracted self-diffusion coef- 
ficients as a function </>. To estimate an ideal glass line, 
we apply an MCT-inspired fit D ~ (<f> — g ) 7 , for the 
two polydispersities at which crystallization never inter- 
venes. The data for <fr = 0.59 are excluded because these 
samples exhibit ageing (we discuss ageing further below). 
The fits give 7 ~ 2.21 and 4> g ~ 0.586 (for s = 0.085) 
and 7 ~ 2.15 and 4> g ~ 0.585 (for s = 0.07). Both fits 
are included in Fig. [2] (bottom) ; their difference is barely 
significant. Using the computed D values, we next es- 
timate isodiffusivity lines in the cj>, s plane. These lines 
can be measured outside the crystallization region, and 
also some way into its interior (in cases where samples 
have exited the MSD plateau before they crystallize). 
The isodiffusivity lines shown in Fig. [1] are parallel as 
for calculations for other repulsive systems [lo, HH , and 
almost vertical (s-independent). We also show in Fig. [1] 
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FIG. 3: MSD versus t - t w for $ = 0.60 and s = 0.06, at 
various t w , averaged over 9 different runs. The MSD for the 
same (j> and lower s (s = and 0.03) is also included. At these 
low s the system crystallizes, and averaging over different runs 
is no longer possible. As a consequence, the MSD is more 
noisy but, within the statistical error, identical to that of 
s = 0.06 until crystallization takes place. The dashed-dotted 
line, with slope 1, has been included to stress the fact that 
the system is subdiffusive at the moment it crystallizes. Inset: 
Pressure versus time for the same state points. Vertical lines 
indicate the i^-values employed in the MSD calculations. 



FIG. 2: (top) Time evolution of the average MSD for 4> = 
0.54, 0.55, 0.56, 0.57, 0.58 (from top left to bottom right) for 
all non-crystallising state points. Most symbols are indistin- 
guishable due to the almost perfect collapse of data for dif- 
ferent s. (bottom) Extracted D values versus (j> (symbols), 
with corresponding power-law fits for s — 0.07 (dashed line) 
and for s — 0.085 (dashed-dotted line). D values at <f> — 0.59 
show ageing (see below) and are excluded from the fit. 

the glass line from an MCT calculation for a binary mix- 
ture with s at each point equal to that of our simulations 
using as input the partial static structure factors 
calculated from additional simulations. On moving away 
from the monodisperse limit ((f>MCT — 0.525 this is 
near-vertical and tracks the isodiffusivity and extrapo- 
lated ideal glass lines. 

The above data presents strong evidence that polydis- 
persity up to 0.08 has barely any effect on the dynamics 
of structural rearrangement in amorphous HS systems 
on the approach to arrest; independently of polydisper- 
sity, these show a glass transition (ideal or otherwise) 
at <f> = 4> g « 0.585. This points firmly toward destabi- 
lization of the crystal as the reason for polydispersity to 
suppress crystallization. Remarkably though, Fig. Q] also 
shows that systems of low polydispersity (s < 0.05) can 
crystallize for <j) well beyond the extrapolated 4> g . 

The most obvious explanation for this would be that 
these systems are after all ergodic: the glass transition 



is not ideal. What matters though is whether they are 
ergodic on the simulation timcscalc. To check this we 
have studied closely the case of <f) — 0.6, for both crystal- 
lizing (s < 0.03) and non-crystallizing (s > 0.03) poly- 
dispersity. In Figure [3] we report the MSD as a func- 
tion of t — t w for <j) = 0.60 and s — 0.06, where t w is 
the ageing time elapsed from the beginning of the run. 
While the short-time and intermediate plateau regimes 
remain essentially unchanged, a significant ageing of the 
dynamics appears at long times with increasing t w . In 
the inset, we report the pressure versus time for the same 
state points. For <p = 0.60 and s = 0.06 the pressure 
does not equilibrate to a clear plateau, but continuously 
decreases: hence this system is nonergodic on all time 
scales simulated (t < tM = 10 5 ). Invoking now our 
previous evidence that polydispersity has no discernible 
effect on dynamics in the amorphous state, we can in- 
fer that the same nonergodicity is present in amorphous 
states at <f) = 0.60 and s < 0.03, which do crystallize. 
Indeed, by monitoring the MSD in the MD simulations, 
we find that these near-monodisperse samples crystallize 
on time scales ~ 10 3 (shown in Fig. 3 by a sharp rise 
of the MSD and a sharp pressure drop), at least two or- 
ders of magnitude smaller than the timescale to attain 
ergodicity (inferred from the ageing studies). Fig. 3 con- 
firms that the MSD curve found (on a fresh sample) at 
high s is closely tracked by the corresponding low-s sam- 
ples, until interrupted by crystallization. Moreover, on 
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these crystallization time scales, particles have diffused, 
on average, only a few tenths of a diameter. 

We conclude that dense hard spheres can crystallize on 
a time scale short compared to that of particle diffusion in 
the initially homogeneous, amorphous state. Put differ- 
ently, bulk hard sphere glasses can crystallize. The mech- 
anism of this transition, direct from amorphous solid to 
crystal, is currently under study. A strong possibility, al- 
ready proposed in [3], is that this direct crystallization 
channel is closed only when the terminal polydispersity 
is reached. Only beyond this point is diffusion needed, to 
allow fractionation to occur. (Suggestively, the crystal- 
lizing region in Fig. [T] lies wholly within the unfraction- 
ated crystal zone on the equilibrium phase diagram of 
jUUHl-) Note also that, to equate pressure, the crystal 
is denser than the glass from which it forms [Io| . 

Our findings undermine an implicit argument in the 
colloid literature @, 0, E3] that if an amorphous state 
readily forms crystals in bulk, this state must have been 
an ergodic fluid, not a glass. This now appears unsustain- 
able, and one can no longer infer [13], from the observed 
rapid crystallization at <f> > cf> g of colloids in microgravity, 
that such colloids have no glass transition. The view that 
crystallizaton implies mobility is anyway rather at odds 
with another longstanding view, that truly monodisperse 
systems crystallize on a timescale so fast that a glass can- 
not be seen (see [ijjj] for a discussion) . This view also ap- 
pears overstated; our results for cf> = 0.63 and s — show 
glassy dynamics without crystallization up to t ~ 10 5 (an 
experimentally accessible range for colloids). 

In this Letter we have studied by simulation the in- 
terplay between dynamics and crystallization in both 
monodisperse and polydisperse hard spheres. We have 
computed a stability diagram in the volume fraction / 
polydispersity plane, showing the region where crystal- 
lization takes place. This region is controlled at low <j) 
by supersaturation (the thermodynamic driving force for 
crystal nucleation) and at high <j) by an interplay be- 
tween polydispersity and dynamical arrest. Computing 
the mean square displacement, we concluded that poly- 
dispersity barely affects the dynamics of diffusive rear- 
rangement in the amorphous state. An MCT extrap- 
olation of the isodiffusivity lines yields, for all studied 
s (including s — 0), an ideal glass transition point at 



ported by eDIKT. 



0.585. For 



> 



the system is nonergodic, and 



shows ageing. However for near-monodisperse samples 
this process is interrupted by a glass-to-crystal transition, 
whose mechanism we hope to elucidate in a future work. 
We conclude, crucially for colloid experiments, that for- 
mation of a crystal does not prove absence of a glass. 
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